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During fabrication and scanning-tunneling-microscope (STM) imaging of dangling 
bonds (DBs) on a hydrogenated siHcon surface, we consistently observed halo-like fea- 
tures around isolated DBs for specific imaging conditions. These surround individual or 
small groups of DBs, have abnormally sharp edges, and cannot be explained by conven- 
tional STM theory. Here we investigate the nature of these features by a comprehensive 
3-dimensional model of elastic and inelastic charge transfer in the vicinity of a DB. Our 
essential finding is that non-equilibrium current through the localized electronic state 
of a DB determines the charging state of the DB. This localized charge distorts the elec- 
tronic bands of the silicon sample, which in turn affects the STM current in that vicinity 
causing the halo effect. The influence of various imaging conditions and characteristics 
of the sample on STM images of DBs is also investigated. 
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Dangling bonds (DBs) on silicon crystal surfaces are becoming increasingly important can- 
didates for building extremely miniaturized nanoelectronics devices. DBs are band gap states 
that can act as donors or acceptors, depending on their energy level with respect to the Fermi 
level. A surface silicon atom can host a dangling bond when only three of its valence electrons 
are involved in covalent bonds, the remaining valence electron being left in an unsatisfied 
(dangling) bond state. It turns out that a DB can be in three, energetically different states: 
a negative state DB^ (holding two localized electrons); a neutral state DB^ (holding a single 
localized electron); and positive state DB^ (holding a localized hole). 

The exact energy levels for an isolated DB on an otherwise perfectly terminated H:Si (001) 
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2x1 surface are not exactly known, but reasonable estimates can be obtained by various com- 
putational methods: an extended Huckel theory lUl, Poisson-Schroedinger equations HI, and 
density functional theoretical (DFT) methodslISj. The three-dimensional form of the DB wave- 
function must, as a surface state, decay exponentially into vacuum and the crystal, as well 
as laterally away from the center of the host Si atom. Theoretical estimates lU-ISl and experi- 
mental measurements (by contactless capacitance-voltage method! 4 1) of the DB-state energy 
in the bandgap are found to agree relatively well. In this study we assume the energy levels 
are, with respect to the silicon valence band maximum (VBM): £'db-= 0.82 eV [|4|, £"060= 0.35 
eV O, £^DB+= 0.0 eV 121. Because the latter state is degenerate with the silicon valence band, 
in this study it is assumed not to play a significant role. 

Studying constant-current empty-state scanning tunneling microscope (STM) images for 
different sample bias voltage, Vbias, and current setpoints reveals in some cases the existence 
of a dark halo surrounding a bright spot of atomic size at the location of the DB. The appear- 
ance of such a halo has been also observed when imaging single dopant atoms located near 
the surface for other types of semiconductor Il6l [71. The simple and intuitive explanation is 
that this dark halo is caused by upward bending of the energy bands at the surface, lead- 
ing to a reduction of the STM current (in constant-height imaging mode) or equivalently a 
decrease in the STM tip height (in constant- current imaging mode). In order to distinguish 
such spatially limited band bending from one of a large scale, e.g. for a great density of DBs, 
it is important to understand to what extent band bending due to the charging of a collective 
of such surface states occurs uniformly along the surface and where the delimiting behavior 
occurs. Arguably, a rule of thumb is that band bending becomes uniform (and the Fermi level 
becomes pinned at the surface) when the average distance between two dangling bonds be- 
come equal to the screening length in the semiconductor. However, for the cases approached 
in this study, due to the great extent of H- termination of our silicon surfaces, this condition 
is not fulfilled. Therefore, the Fermi level is not pinned at such surfaces. This fact implies 
that each isolated DB (unless very close to others), has its own individual screening field, 
due to mobile charge carriers in the crystal. On the other hand, isolated groups of closely 
spaced DBs (< 2nm) are electronically coupled IIH |3l |8j, and have more or less a common 
screening field around the whole DB group. The coupling of such groups in the presence of an 
STM tip was recently studied experimentally and theoretically within an extended Hubbard 
model IIHJ to yield time-average occupations of individual DBs, without explicitly accounting 
for the dynamics of electrons in/out the DB group. 



(a) 



(2V, 1nA) 




(2V, 0.2nA) 



(2V, 0.1 nA) 



FIG. 1: STM images of a 4.0 x 3.6 nm region of the 2x1 H-Si (001) surface with an isolated DB in the 
unoccupied imaging mode (sample biased positively) at imaging setpoints (sample bias and current) 
indicated below each image. 



The properties of various types of single defect levels in semiconductors have been pre- 
viously studied in the context of their STM imaging appearance |[9|. For DBs, the carrier 
capture properties llTOl [Till and the on-site Coulomb interaction betvi^een two localized elec- 
trons can be extracted from spectroscopic measurements [il2i] on isolated DBs on a B-doped 
Si(lll)-(A/3 X y/3)R30°. However in these latter studies, electron dynamics is simpler than 
for a DB on hydrogenated Si surface, as there is no tip-induced band bending at the surface. 
Here we study single electron dynamics for the case when band bending is present at the sur- 
face. We account for all the important charge transfer mechanisms in order to gain insight 
into the charging dynamics of a DB and its appearance in STM images. 



II. DBS IN THE POTENTIAL LANDSCAPE OF THE IMAGING PROBE 

In Figure [l] we show typical unoccupied-state STM images of dangling bonds on H- 
terminated Si (001) with a 2x1 surface reconstruction. In part (a) both single and coupled 
DBs are seen, with halo-like features. The halos around single DBs are consistently more 
prominent than those around groups of coupled DBs. In (b) a single DB can be identified at 
the center of the image as a bright spot surrounded by a darker ring or halo-like feature. This 
appearance is consistent with many experimental observations taken at similar conditions, 
and has been reported in the literature lUHl. Also, such imaging features have been observed 
in the case of other localized charge centers, such as subsurface dopants |]9j. 
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We define an isolated DB on the silicon surface as one with no equivalent nearby DBs to 
be significantly coupled with (nor other resonant electronic levels). This implies that at an 
isolated DB an electron is strictly localized in an orbital at a single Si atom. We assume a 
simple Slater-type orbital (STO) for the DB wavefunction having p-type symmetry (two lobes) 

V'DB(r) = ^DB{r, e) = Nr cos{e) exp[-C(0)r], (1) 

where r and 6 are the three-dimensional spherical coordinates and is a normalization con- 
stant. The difference from a "regular" STO is that in our case the decay rate C in the exponen- 
tial must vary with the z-coordinate in order to ensure consistency with the asymptotic decay 
in the specific environment of the DB, as follows. The decay rate is related to the ionization 
potential Wi of the DB electron via C = \/2rnWJh. This ensures the correct asymptotic be- 
havior (from first-principles) of the wavefunction above. As the DB orbital is located partially 
in vacuum and partially in silicon, the ionization potential is not a constant, but it rather 
depends on the "ionization path" of the DB electron, e.g. it is energetically easier to extract 
an electron toward the bulk (into the conduction band of silicon) than toward the vacuum. To 
reflect this fact, here we assume a simple form for this dependence, namely 

W,{z) = lybuik + \ (H^vac - VFbulk) [tanh(z/u;) + 1] (2) 

where Wbuik and W^a.c are the ionization potentials with respect to bulk and vacuum, respec- 
tively, z is measured from the surface and positive toward vacuum, and w is a characteristic 
width of the transition as can be seen in Fig. [2] (a). The normalization of this wavefunction 
is calculated numerically and the resulting modified STO function is plotted in Fig.[2](b) as a 
color map for a two-dimensional axial section. Note that, as a consequence of our choice for 
the C function, the two lobes of the DB orbital are quite disproportionate, with the vacuum 
lobe being much smaller in spatial extent and portion of charge. This is also qualitatively 
consistent with DB orbital characteristics derived from ah initio calculations for an isolated 
DB at the surface of model hydrogenated silicon clusters[3|, mimicking the 2x1 H-Si (001) 
surface. 

When the DB is negatively charged, the localized extra electron is expected to locally dis- 
tort the electronic levels at the surface of the host silicon crystal, i.e. cause band bending. 
Band bending in a semiconductor is the phenomenon of distortion of electronic levels of the 
crystal in the presence of an external perturbing electric field. During STM imaging, another 
source of band bending is the biased STM probe itself and this is known as tip-induced band 
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FIG. 2: Characteristics of the DB orbital, (a) The ionization potential as a function of distance from 
the surface (positive z in the vacuum), (b) Color map of the distribution of charge density of a DB in 
an axial plane as predicted by the modified STO function in eq. [l] The color bar indicates the values 
of the density in units of 10^'' m^'^. 

bending. Depending on the value of the STM bias voltage, the amount of band bending can be 
considerable, sometimes enough to induce an inversion layer near the surface (whenever the 
Fermi level becomes lower than the intrinsic level of the semiconductor). The amount of tip- 
induced band bending can be theoretically estimated by solving the Poisson and Schroedinger 
equations simultaneously, for example by finite element methods (FEM) for an assumed probe 
geometry. 

During STM imaging of isolated DBs, as the imaging tip is brought closer to the DB loca- 
tion, the band bending is turned gradually up and so is the DB level relative to the sample 
Fermi level. In Figure |3j we show a typical band bending diagram calculated by FEM for the 
silicon surface in the presence of a scanning probe tip. The sample is n-type with a resistivity 
of 3.5 mJl-cm, the sample bias voltage was assumed +2 V, and the tip height was assumed to 
be 7 A. For these typical imaging conditions, the bands are bent upward to the extent that the 
DB level is higher than the sample Fermi level, by 0.3 eV. This value is much greater than 
kT, which - assuming thermodynamic equilibrium between the DB state and the bulk semi- 
conductor - means that at room temperature, the DB state should be completely unoccupied. 
However, the process of STM imaging also inject electrons into the DB state and an interplay 
of in/out transfer rates starts to form. In this study we show how this fact has important 
consequences on the way surface states are imaged by STM. 

An important consequence is evident from the diagram in Figure [3| in the presence of the 
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FIG. 3: Calculated band bending diagram as a function of the z-coordinate. The black curve is the 
CBM, the red dotted line is the sample Fermi level, and the blue line shows the DB^ level with 
indication to tunneling into resonant levels in the silicon conduction band. The sample is n-type with 
a resistivity of 3.5 mfi-cm, sample bias +2 V, STM tip height 7 A. 

tip-induced field, the DB level becomes resonant with CB levels of the host crystal, implying 
that tunneling from DB into the crystal becomes possible. This situation is then similar 
to double barrier tunneling junctions, which were recently used to study energetic levels of 
single atoms Ifl4l . defects IflSll . and molecules lfT6] - [l8]| . at thin insulating films on metallic 
surfaces, as well as metallic nano-islands on Si lHOll . Manipulation of charge states llT6ll and 
bond structure Il20ll of molecules has also been demonstrated. However, the similarity of our 
system to these other cases is limited because in our case there are a number of other factors 
at play that contribute to the charging of an individual DB, as we show below. 

Careful examination of our available experimental data in conjunction with a quantita- 
tive theoretical analysis, indicates that STM imaging features of dangling bonds are due to 
dangling bond charging caused by non-equilibrium effects during imaging, i.e. rate-limited 
charging/discharging of a DB during extraction/injection of electrons by the STM tip. For 
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routine unoccupied-state imaging conditions (sample bias +2 V, current setpoint lOOpA- InA) 
a DB state in electrochemical equilibrium with the bulk silicon would be unoccupied un- 
der the STM tip, but because the rates of charge transfer are limited, the DB state fails to 
achieve equilibrium and becomes instead a non-equilibrium steady state. STM electrons in- 
jected into the DB become temporarily trapped there during imaging due to a slow discharge 
mechanism into the Si sample, which cannot keep up with the rate of injection by the STM 
tip. Explaining the DB imaging as purely quantum mechanical behavior, namely as Friedel 
oscillations in the local density of states has been proposed - without consideration of band 
bending and nonequilibrium effects - but this explanation is inconsistent with the variation 
of the DB features with temperature and STM bias voltage II2TII . 

III. NON-EQUILIBRIUM CURRENT THROUGH A DB STATE 

In the above discussion we assert that a net negative charge on the DB during imaging 
can only be explained by a non-equilibrium phenomenon, and is in fact impossible if the DB 
was in electrochemical equilibrium with the silicon crystal (due to the DB level being greatly 
elevated from the bulk Fermi level). In our system, we have a situation where the current 
flows from one reservoir (metal) at the chemical potential /itip = fJ-i to another (Si crystal) at a 
chemical potential /isi = ^2 via a discrete level (DB) whose energy is somewhere between /itip 
and ^si- An electron can be transferred into or out of the localized DB level at two different 
rates, say rtip = ri and rsi = r2, which in general include all elastic and inelastic charge 
transfer mechanisms. Overall, a steady state of current flow is achieved when rtip = ?^si in 
which the occupation of the DB level depends on the DB energy and the escape rates. 

For a DB state in equilibrium with reservoir 1 only (or 2 only), its occupation equals 
the Fermi-Dirac distribution, f^l = /fd(^db,^i) = [1 + exp((£^DB - ni)/kT)]-'^ (or fjyl = 
/fd(-Edb, = [1 + exp((£'DB — fJ'2) / kT)]^^ , respectively). However, since the DB level is in 
contact with both reservoirs, the DB will achieve a steady state of filling, with an occupa- 
tion /pg somewhere in between (the superscript * stands here for non-equilibrium occupa- 
tion function). Then, the current through one contact equals the current through the other, 
Ii = h, which is the condition for a steady state occupation of a DB. In order to calculate the 
steady-state current and DB occupation, we must first determine the flow rates between the 
tip and DB on one hand, and between sample and DB on the other. The possible mechanisms 
of charging and discharging a DB in the empty state STM imaging mode are shown in Fig. |4} 
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FIG. 4: Three-dimensional rendering of an isolated DB on the H-Si(OOl) surface in the presence of 
an STM tip in the unoccupied imaging mode, (a) Various channels for elastic and inelastic charge 
transfer between the STM tip and the silicon sample in the vicinity of the DB are indicated by colored 
arrows and labeled for further discussion. Si atoms are green, H atoms are white and the DB is shown 
schematically as purple sphere, (b) Diagram of the silicon band gap, which further explains the various 
channels shown in (a): (2) inelastic capture by the DB° level (black horizontal line) of an STM electron 
after which DB becomes negative; (5) thermal electron emission from the DB~ level (red horizontal 
line) into the CB; and (6) inelastic capture of a hole by the negative DB, followed by recombination 
and the DB becoming neutral, (c) Tip-induced band bending diagram in the direction normal to the 
surface. The elastic charge transfers shown in (a) are detailed here: (1) direct tip-silicon tunneling; (3) 
tip-DB tunneling across the vacuum gap; and (4) tunneling from the DB~ state into the CB across the 
potential barrier created by band bending. The Fermi level in silicon is indicated by a dashed red line. 
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Note that, although in principle the DB+/DB° transition should also contribute to the 
transfer of electrons from the tip to the silicon crystal, its contribution can in practice be 
neglected for subtle reasons derived from the fact that the DB+ level is resonant with the VB 
band (and therefore the DB+ state has a negligible lifetime). As a consequence: (i) tunneling 
of electrons from the STM tip to the DB^ level almost never happens, because it implies 
a preexisting DB+ state; (ii) thermal emission of an electron from the DB'^ level involves 
an additional thermal barrier of 0.45 eV (or IS/cb^), compared to emission from the DB~ 
level, and it can be completely neglected. Therefore, the function /^g signifies the statistical 
occupation of the DB^ level alone, and varies between (for a DB with a time-average of zero 
charge) and 1 (for a completely filled DB , on average). In the sections below we discuss each 
charge transfer mechanism shown in Fig.|4]in detail. 



A. Rate of electron transfer between STM tip and an isolated DB 

During a STM experiment, a possible mechanism for electron transfer from the STM tip 
to a DB (and vice versa) is by elastic tunneling across the vacuum gap. Another possible 
mechanism is a two-step capture process in which an electron from the tip first tunnels into 
the Si crystal in the vicinity of the DB and then immediately gets captured by the DB before 
it has the chance to become an extended- state electron. Henceforth we call this mechanism 
"STM vicinity capture of tunneling electrons by a DB". 

1. The rate of elastic transfer between the STM tip and an isolated DB can be estimated 
(see Appendix |A] for more details) 

4:7r€ 1 ~ 

-^tip-DB(d) = — Ptip(-EDB)[/tip(-E'DB) — /DB]|-^tip-DB(-E'DB, d)| , (3) 

tl i^DB 

where we denote the Fermi-Dirac distribution simply by / with a subscript indicating 
the respective electrode ("tip" for tip, or "Si" for silicon crystal), JIdb is the volume of 
the DB (defined as the region in space where the DB charge density is greater than 1% 
of its maximum), and |M(£'db, d)P is a properly normalized transfer matrix element for 
separation d (between DB and tip centers) defined as 



|Mtip„DB(^DB,o!)P 



2 \ 2 
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(4) 



Note that in eq. [3] above, we also replaced the Fermi-Dirac distribution function for 
the DB with the non-equilibrium occupation function. As explained in Appendix |A[ the 
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functions above represent dimensionless wavefunctions of the respective states under 
consideration defined as 

ifir) = VUi^ir), (5) 

where V is the normahzed wavefunction and Q is the volume of the electrode or the 
localized state, respectively. For the purpose of estimating the tunneling currents, the 
wavefunction of the STM tip is taken in the conventional form of the Tersoff-Hamann 
approach 

for |r — rtipl > Rq, where rtip is the position of the tip center, Rq is a characteristic 
tip radius, ctip(i?o) = cotCtip-^o exp(Ctip-Ro)/\A^j cot is a dimensionless constant of the 
order of unity, Ctip(-S) = \j2m{E^^c — E)/h, and E'vac the vacuum level. 

Generally speaking, one would expect sequential tunneling events from the STM tip, 
first into the DB° level and then into DB level. However, for an electron to tunnel 
into the DB^ level, the DB has to be in a pre-existing DB+ state, and the probability of 
having that state is negligible because of it being resonant with VB levels. Therefore 
tunneling into the DB'^ level can be neglected. 

2. The STM vicinity capture of tunneling electrons by a DB involves tunneling from the 
STM tip into the Si crystal for a range of electron energies and momentum orientations 
and subsequent capture of those hot electrons by a neutral DB. The corresponding cur- 
rent into the DB is denoted by C'^^ and the details of its calculation are described in 
Appendix [Cj 

During an STM experiment, in addition to the tunneling current from the tip to the DB, 
other transfer mechanisms contribute to the overall balance of charge on a DB orbital. They 
are discussed in detail below below. 

B. Rate of electron transfer between an isolated DB and the silicon crystal 

Once an incoming electron becomes localized at a DB state, there are three main mecha- 
nisms of transfer into the silicon crystal: 
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1. The localized DB electron escapes by emission into the conduction band via thermal ex- 
citation. This thermal escape rate only depends on the barrier height and temperature; 
the escape rate is given by 

fin = CJnfnA^cexp[-(£^CBM " EBB)/kBT], (7) 

where the prefactor is an attempt frequency and Nc is the effective density of states at 
the bottom of the conduction band 

1 f m^cskBT'^^^'^ 



with mn,efr being the effective mass of electrons in Si. 

2. The localized DB electron tunnels elastically into the bulk CB level (at the classical 
turning point level where the CBM intersects the DB level). Henceforth, CBM stands 
for CB minimum, and VBM stands for VB maximum. Tunneling from DB occurs across 
the space charge layer and the tunneling distance depends, via band banding, on the 
doping level of the sample, STM tip height and bias voltage. The tunneling current can 
be estimated by using the formula (see Appendix^ 

AlTF 1 

/DB-Si = ^^5Si(i^DB)[/5B - /si(i5^DB)]|MDB-Si(i^DB)P, O) 

where Mdb-si in this case also depends on the tunneling barrier for an electron with 
energy Edb which is given by the band bending dependence on z. Upon tunneling into 
the bulk, the DB resumes a neutral state. 

To complete the picture, for the purpose of calculating transfer coefficient between the 
tip and the silicon sample, we assume sample wavefunctions derived from the "jellium 
model"! 22] with modifications to account for the corrugation of the surface 



V'f kW = -csi ^ ; — exp(iK • R), (10) 

where E is the eigenenergy of the state; K = {kx,ky) is the surface parallel wavevec- 
tor; kz is the surface normal wavevector (the last three being related via E = 1t?{K'^ + 
kl)/2m); r = {x,y,z), R = (x,y); i7si is the sample volume; Cli = 2m(£^vac — E)/h?; 
and csi = icos/\f^i with cqs a dimensionless constant of the order of unity. The corru- 
gation function of the sample H(R) is in general two-dimensional, but in our model 
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we treat the x and y directions as independent. We assume the corrugation com- 
ponent along dimer rows (for which we calculate STM current traces) has the form 
H{x) = Acos^{'Kx/Lx), with A the corrugation amplitude and its period. 

3. A hole striking the surface is captured by the occupied DB level (defect mediated re- 
combination). The associated recombination current is calculated in a customary ap- 
proximation, by assuming that the rate of recombination at a DB site is proportional to 
concentration of holes, p, at the vicinity of the DB, the thermal velocity of those holes, 
Vp, and a capture cross section, of a hole by a negatively charged DB, a^: 

Tree = CTpPVp (11) 

Their numerical values are Vp = 1.87x 10^ cm/s, dp = 1.1 A^, while p depends on the band 
bending at the surface and can be anywhere from 70 cm in the bulk of highly doped 
n-type Si (resistivity 0.01 il-cm) to 10^^ cm in an inversion layer. The recombination 
current at the DB is given by 

/rec = Cp = eCTpPfp/oB (12) 

An additional note of the capture cross section: defects in semiconductors, both in bulk and 
at surface/interface, have been known to act as very effective recombination centers. In fact, 
carrier recombination in a semiconductor is actually dominated by trapping processes at de- 
fects, such as DB states. An important fact is that defects with energies close to band edges 
(as opposed to the middle of the band gap) are more efficient at trapping carriers and causing 
recombination [|23i . Also, the above rough values of recombination rates rely on the semiclas- 
sical estimates for the hole concentration at the surface and include no charge quantization 
effects due to band bending. This quantum effect can significantly alter the carrier concen- 
tration at the surface and therefore can be important for the recombination rates. 
Generally speaking, the capture cross section has a thermally activated behavior, 

<7p = f7opexp(£^a/A;Br) (13) 

where aop is a constant and Eg, is an activation energy. For our purpose, we note that the 
capture cross section of a DB is more or less independent of temperature between 100-400 K 
lEH, for both electron and hole capture. It follows, that in the absence of tunneling, the flow 
of electrons in/out of a DB is given by 

Rn = Cn- En = eanVnn{l - /db) - eCn/DB (14) 
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with capture competing with emission, and the unknown occupation, /^g dictated by the 
non-equilibrium conditions. 

IV. SOLVING FOR THE NON-EQUILIBRIUM OCCUPATION OF A DB 

Combining all the above mechanisms we obtain the total current flow between DB and 
bulk 

-^DB = -^D^) (15) 

yields 

-^tip-DB + Cn™ + Cn = -^DB-Si + -E'n + C'p, (16) 

where the additional current components, in Ampere, are given by 

Cn = eC7nfnn(l - /db), (17) 

En = een/^B, (18) 
C^™ = ea„/^>B_°g,(l - /Sb) / dr|^DB(r)|^^^, (19) 
Cp = eapVppf^ji, (20) 

where Vp is the thermal velocity of holes and Ny is the effective density of states at the top of 
the valence band. 

_ 1 fmp^skBTV^^ 1 f mhp.efffcsr y/^ 

with mip eff and mhp eff being the effective masses of light and heavy holes in Si. 
The solution is 

itip /tip(^DB) + ^Si/si(^DB) + Cn + 4™ 
Hip + isi + Cp + Cn + Cn + d 

where the charge transfer rates, in s~^, are given by 
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(23) 
(25) 

<7tip(^DB)|Mtip-DB(£^DB)|', (26) 
5Si(£^DB)|MDB-Si(£^DB)|', (27) 
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and where D{9{r)) in c„ is the normahzed tunnehng current density in the direction 9 from 
the STM tip to the sample, which is derived in Appendix C. 
The total time-average STM current is then given by 

^STM = ^i?un?(l - /db) + -^SST/dB + -^DB: (28) 

where I^^^ and I^^]^ are the tip to sample currents in the presence of a neutral and negative 
DB, respectively. 

From ab initio calculations on isolated DBs at the surface of model silicon clustersfSl, we 
know that a negatively charged DB renders its host Si atom to be slightly elevated from the 
plane of the surface, about 0.3 A, which can also have an effect on the tunneling coefficient 
during STM imaging. This effect is also captured in our model by accordingly modifying the 
sample wavefunction, and plays a role in the appearance of the DB. 

Another possible effect during DB imaging is a Stark shift of the DB level due to the tip- 
induced field in which the DB is found. This shift would be most pronounced when the field 
is stronger, i.e. when the tip apex is immediately above the DB, for a given sample bias. In 
principle, calculating the magnitude of this shift is possible within the frame of our current 
theory by solving the Schroedinger equation in the presence of the external field, e.g. by 
a perturbative approach or finite element method. This would produce another parameter 
in our model, namely the magnitude of the Stark shift as a function of the tip-DB distance 
and tip height. However, for the sake of simplicity, we decide not to include this additional 
parameter in our model, but rather to comment on its possible consequences on the model 
results. 

Note that if the "equilibrium would-be" populations are equal the net current through the 
DB is zero (this only happens if the chemical potentials on the two sides are equal). Also note 
from above formulas that, if one of the flow rates is much smaller than the other, then the net 
current will be reduced to just that rate. As another general remark, an implicit assumption 
above is that the total STM current (through DB and bulk states) is small enough that it does 
not significantly change the Fermi level of the reservoirs in the vicinity of the two contacts. 

V. RESULTS OF THE IMAGING MODEL 

In this section we present the numerical results of the theory presented above for the 
imaging of DBs. We illustrate our model with experimental and theoretical results on a 
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highly doped n-type silicon sample, with a H-terminated 2x1 (001) surface reconstruction. 
We only focus on the unoccupied-state imaging mode of STM, where most of our experimental 
data was obtained. In this case, upward band bending occurs at the surface and an excess 
concentration of holes is present there. The STM tip injects electrons into the sample, and into 
isolated DBs as well. Localized electrons at DBs discharge into the sample by the different 
channels described above, including recombination with the excess holes at the surface. Other 
sample doping type and levels are also amenable to this model, keeping in mind that the 
system parameters will be different and require making appropriate adjustments: Fermi 
level, band bending, local electron and hole concentrations at the surface, etc. For example, 
in the occupied-state imaging, electrons are extracted from the sample into the STM tip, and 
downward band bending occurs at the surface, which brings about an excess of electrons, as 
opposed to an excess of holes. This has consequences on the filling of isolated DBs, prior to 
their discharge to the STM tip. 

Finite element analysis was used for calculating tip-induced band-bending as the self- 
consistent solution of Poisson and Schroedinger equations Il25]| . This solution is reliable for 
the regime in which quantum confinement effects (quantization of electronic levels at surface) 
are negligible, but also for the case of strong inversion or accumulation, when quantum effects 
are import ant ll26ll . In this latter case, many surface quantized levels are occupied, which 
makes the semiclassical approximation for the charge distribution valid. The semiclassical 
calculation fails to be reliable when the inversion or accumulation is moderate, that is, when 
surface quantized states exist but only a few quantum states are occupied. Some details of the 
FEM calculation are given in Appendix |B] and the complete details of the whole calculation 
procedure are described elsewhere I125II . 

In Figure [5] we show an example of calculated potential landscape (by FEM) for chosen 
sample bias and tip height as (a,b) a two-dimensional section of the axially symmetric system, 
and (c) a function of the depth inside the silicon surface. The effect of a negative DB under the 
tip on the band bending was also calculated and shown here. The greater differences between 
the two cases can be seen in the silicon region immediately under the tip apex, where the DB 
is located. Simultaneously, we calculate the carrier concentrations in the sample induced by 
the tip alone and also by the tip in the presence of a DB and plot the results in panels (d,e). 

Unlike the case of a highly doped p-type Si sample, where the main discharging rate of a 
DB~ is that of recombination with holes IITOl . in our case, both types of carrier concentration 
are small at the surface, where the sample is rendered almost intrinsic by the tip-induced 
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X [nm] 

FIG. 5: (a) Two-dimensional section of the electrostatic potential landscape caused by a biased STM 
probe in the vicinity of a silicon surface for unoccupied state imaging and for axial symmetry (tip 
central axis is at a; = 0). The potential penetrates the surface and vanishes deep in the crystal. The 
Si crystal is n-type with a resistivity of 3.5 mQ-cm, voltage bias on the STM tip is -2V, and the tip 
height was assumed to be 0.7nm (tip apex is at z = 0). The tip contour is drawn as a black line and 
the surface as a white line, (b) The electrostatic potential for the same biasing conditions as in (a) and 
in the presence of a negative DB directly under the tip apex. The color bar indicates the values of the 
electrostatic potential in V, and is identical for both plots for easy comparison. (c,d,e) Subsurface band 
bending and carrier concentration along the central axis of the tip plotted as a function of sample 
depth, z, for biased tip in the presence of a DB° (black curves), and in the presence of a DB (red 
curves) placed directly under tip apex. 

potential. As a consequence, recombination is no longer the dominant discharging rate and it 
competes vi^ith the other rates described above. 

The calculated band bending and carrier densities are input to the mathematical formal- 
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ism described above in order to compare our theory to our experimental observations. As 
suggested above, the key point is to determine the charging state of the DB. The main results 
needed to understand the non-equilibrium charging effect are shown in Figures [6] and [7| Part 

(a) illustrates the competitions between electron transfer rates in and out of an isolated DB in 
the unoccupied-state imaging mode. As the STM tip moves in from far away toward the DB, 
the electron injection rate from tip to DB (solid black curve) increases approximately as an 
exponential function, up to a critical distance Tc of about 5 A. At that point, this injection rate 
becomes greater than the dominant rate of transfer out of DB, which in this case is provided 
by the thermal emission of electrons into the conduction band. The second greatest discharge 
rate is by tunneling into resonant CB levels, and the third is by capture and recombination 
with a itinerant hole, Cp, which is considerably smaller. 

For smaller distances than rc, the discharge rate becomes the limiting rate for current 
through the DB, and therefore it dictates the charging state of the DB, /j^g, plotted in part (b). 
In turn this charging state determines the total band bending (tip-induced plus DB-induced) 
and also the tip-sample tunneling coefficient in the vicinity of the DB and it therefore affects 
the total current measured in the STM scan, see part (c). Over a relatively short distance, 
in the vicinity of r^, the time-average DB occupation switches from zero to one, because the 
electrons localized there cannot escape at the same rate at which they are injected. With the 
addition of a localized negative charge at the DB center, the electrostatic landscape changes 
in its vicinity, i.e. upward band bending increases, and the STM current is affected, giving 
rise to a dark halo with the outer edge around rc. this corresponds to the situation depicted 
in the experimental image in Figure [ifa). The bright spot at the DB center is caused by the 
slight elevation of the host silicon atom at the location as mentioned above, which although 
small around 0.3 A, has a significant effect on the tunneling coefficient due to the exponential 
sensitivity of the latter with tip-sample separation. 

A halo does not always appear around an isolated DB, and in our case this happens at 
lower injection current, as in the situation shown in Figure [TJc), which is consistent with a 
greater tip height. To simulate the situation in our model, who performed the same calcula- 
tions for a tip height of 8.5 A and the results are shown in Figure [7| We can see a decrease 
in this case of the STM injection rate into the DB, which becomes lower than the thermal 
emission rate as seen in (a). Therefore the DB is mostly neutral in this situation, as seen in 

(b) . No halo is formed here because additional band bending by the DB does not occur. 
Over the course of many experiments, we also observed other related imaging trends of 
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FIG. 6: Dependence of relevant physical quantities on the surface-parallel distance between the STM 
tip apex and the center of an isolated DB. Voltage bias on the sample is +2V and the tip height is 
7 A from the surface, (a) The dependence of various elastic and inelastic rates of charge transfer 
through a DB in the presence of an STM imaging tip. (b) Non-equilibrium occupation function of the 
DB^ state versus distance, (c) Total STM current dependence on the same distance. This situation 
corresponds to the experimental Figure [ija), where a halo is formed. 

DBs: (i) there is a tip dependence in the exact shape of a DB halo, see Fig. [8| (ii) halo-like 
features around groups of two or more tunnel-coupled DBs. The dependence in (i) can be 
qualitatively described by the fact that different STM tips have atomic size protrusions at 
different locations in the vicinity of the apex. During scanning, as the tip approaches a DB, 
there are different onset points for complete charging (/j^g = 1) as the smallest distance 
between tip and DB (and the corresponding tunneling rate) depends on the location of the 
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FIG. 7: The same quantities are plotted as in figure [6} but for a different tip height, 8.5 A, which 
corresponds to the case in the experimental Figure [ijc), where halo is not formed. 



protrusions for each tip. 

For (ii), the same imaging mechanism is at work for groups of coupled DBs, except the 
locahzed charge is shared among more than one DB, and the charge per atom is less (Si. 
Consequently, the DB-induced band bending tends to be lower Also, due to the shape of the 
STM tip, the tip-induced band bending amounts to different values at different DBs, and 
therefore their levels are pulled upward with respect to the bulk VBM by different amounts. 
The localized charge tends to reside at the state whose energy level is lower (farthest from 
the tip). 

Other possible effects were also taken into account in our analysis. The calculation of 
the charge carrier densities at the surface accounting for charge quantization in the surface- 
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FIG. 8: Dependence of the shape of a DB halo on the particular shape of the STM tip used for imaging. 
The same DB was imaged in (a) and (b), at the same imaging setpoint (sample voltage 2 V, imaging 
current 0.2 nA), but with a tip change between (a) and (b). 

perpendicular direction was also performed in a one-dimensional model (z-dependence only) 
by solving the Schroedinger equation for the quantum well formed by the band bending in 
the z-direction. However, for the situation depicted in Figure |5} band bending is not strong 
enough to induce significant quantization of the space charge layer, and the semiclassical 
treatment is reliable. 

As we mentioned above, during a STM scan when the tip apex is directly above the DB, 
the electric field under the apex can be very strong (around 2 V/nm in the vacuum gap) and 
this can induce a Stark shift of the DB level and a deformation of the DB orbital. This shift 
is of the order of a few tens of meV (from our preliminary estimates) and is not included as a 
parameter in our model. However, this shift has possible consequences on the appearance of 
DB, especially in the region close to the DB center. An upward shift in £^db would decrease 
the band bending barrier and thus increase the isi rate, and also the thermal excitation rate 
Cn (these rates are very sensitive to the exact DB level). This means that more current would 
pass through the DB state itself, contributing to a brighter peak at the DB center We expect 
this effect to become significant only for images taken at lower tip heights (higher setpoint 
currents). 

Band bending in a semiconductor is usually explained by appealing to semiclassical argu- 
ments. However, for the semiclassical theory to be accurate, one must involve the approxima- 
tion that the perturbation in the crystal varies slowly on the spatial scale of the crystal lattice 
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constant. In order to test this approximation for our case, here we also undertake a more gen- 
eral approach that can be used beyond the semiclassical approximation, which can also cast 
insight into the long range behavior of screening of a negative DB. For this purpose, the calcu- 
lation of the electrostatic potential landscape was done by separating the whole problem into 
two components, considered to be approximately independent, according to the source of elec- 
trostatic perturbation: (i) the tip-induced band bending; and (ii) band bending by an isolated 
charged DB alone. Solving these two problems requires significantly different approaches 
caused by the different spatial scales of the problems. A perturbation approach If27] - l29l was 
used for the problem of DB-induced band bending. This approach includes the semiclassical 
approximation and the random phase approximation (RPA) as particular limits, and allows 
us to analyze other quantum mechanical effects, such as Friedel oscillations. 

In order to make the computational task feasible for long spatial range (100 nm), the angu- 
lar dependence of the DB orbital was integrated out in this case (with only radial dependence 
remaining). The DB is considered to be placed in a background charge density determined 
by the tip-induced band bending evaluated by FEM calculations. In Figure [9| we show the 
solution of the (separate) electrostatic problem of screening of a DB by the mobile carriers in 
the semiconductor. The maximum band bending due to the DB , is around 0.17 eV, which is 
slightly larger than the all-semiclassical results presented above. Thermally damped Friedel 
oscillations are indeed present in our results (inset in (a)), but they cannot be associated with 
the appearance of dangling bonds, because: (i) they are extremely small in magnitude com- 
pared to the other effects discussed in our model (by more than four orders of magnitude); 
and (ii) they are manifest at greater distances (> 20 nm) than the length scales associated 
with imaging halos. 

Over the course of many experiments, we noticed a significant effect of the sample prepa- 
ration methods on the way DBs image. In particular, there is a strong correlation between 
the value of the temperature used for flash-heating the sample and the formation of (or lack 
of) dark halos around DBs. This effect is caused by the fact that at high flash heating temper- 
atures (> 1250° C), there is a substantial migration of dopant atoms away from the surface 
region. This dopant depletion happens up to a depth ranging from hundreds of nm to a few 
microns and alters the local Fermi level at the surface. In turn, this has consequences on the 
charge transfer rates discussed above and ultimately on the charging state of the DB under 
imaging conditions. These facts also suggest possible control mechanisms of the charging 
state of a DB, by adjusting the local Fermi level at the surface, e.g. by applying an external 
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FIG. 9: (a) The screened electrostatic potential solution for a negatively charged single DB as a func- 
tion of distance from the DB center calculated using the perturbation approach ||27] - |29l . Inset: the 
long range behavior of the potential displaying thermally damped Friedel oscillations, albeit with 
much lower amplitudes than the potential at close range (the horizontal axis is the distance in nm for 
all plots), (b) Charge density of the mobile carriers (holes) brought about to screen the DB charge. 



potential via a local electrode. This w^as also illustrated and discussed in a recent study on 
groups of closely coupled DBs on hydrogen terminated silicon surfaces [(81. 



VI. SUMMARY AND OUTLOOK 



In this paper we propose a mechanism for understanding the STM imaging of isolated 
DBs on silicon surface, w^ith quantitative focus on the n-type surfaces and unoccupied-state 
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imaging mode. A quantitative analysis of the electrostatics of our system shows that the 
tip-induced band bending plays an important role in determining the charging state of the 
isolated DB by effectively pulling the DB state out of equilibrium with its host crystal. Under 
these conditions, the equilibrium Fermi-Dirac occupation of a DB is replaced by an occupation 
function dictated by the steady-state flow of current through the DB, which is very sensitive 
to the lateral DB-tip distance. Therefore the charging state of the DB can vary quickly upon 
lateral approach of the scanning tip, thus causing halo-like features with sharp edges around 
the DB. The sharpness of the edge is related to the slope of the variation of the DB occupation 
function /^g with the lateral distance, see Figure [6] 

We also found that the all-semi classical approach to band bending by an isolated DB com- 
pares reasonably well, in terms of excess band bending by a DB , to more general approaches, 
specifically the RPA, and that Friedel oscillations do not account for imaging features present 
in this system. 

As the dominant DB discharging mechanism in our specific case study is the thermal emis- 
sion of electrons into the silicon CB, lowering the system temperature should have a drastic 
effect on the DB image, and the halo should be increased. The similarity of our system to dou- 
ble barrier tunneling junction can also be exploited in certain regimes, e.g. by lowering the 
temperature and decreasing the thermal emission, thus rendering the DB-Si tunneling as the 
dominant discharging rate. These will be the subject of a further experimental investigation. 

Our model calculates the total STM current in a time-average view of the system. If the 
time resolution of the STM was high enough, it should be possible to discern a "telegraph 
signal" of the current inside the halo region of an isolated silicon DB, as the band bending 
fluctuates with the charging/discharging of the DB. Our model does not reproduce the STM 
tunneling currents in a strictly quantitative fashion. However it does capture all the relevant 
mechanisms at play in image formation, and the main qualitative trends that were observed 
experimentally, with the variation of system parameters such as doping level, sample bias, 
and current setpoint. Because the charging and discharging rates have exponential sensi- 
tivities to band bending values and distances between electrodes and localized states, their 
particular balance can be greatly affected by small changes in these parameters. Therefore 
imaging of DBs will be different for different semiconductor materials (Ge, GaAs), different 
surface reconstructions (e.g. Si(lll)), for samples doped at lower or higher levels, or p-type, 
and for different tips. The latter is the cause of daily variations in the appearance of DBs in 
the lab. However the general formalism presented here can be adapted to suit those similar 
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A. FORMULAS FOR CALCULATING TUNNELING CURRENTS 

Throughout this study we use a generalized formula ll30l for tunnehng currents between 
two electrodes (1 and 2) of volumes ili and $^2 in the presence of a externally applied voltage, 
VhisLs, equal to the difference between the chemical potentials on the two sides /ii — ^2 

h^2 = ^1 dE[f{E,^i^)- f{E,^i2)]gi{E)g2{E)\M{E)\\ (29) 

Here gi{E) and g2{E) are appropriately generalized local density of states (LDOS) per unit 
energy and unit volume (in units of J^^m^'^) of the two electrodes; \M{E)\'^ is an appropriately 
defined (see below) transition matrix element between states in each electrode having energy 
E. We assume that our system fulfills the conditions in which the above formula is equivalent 
to the golden- rule expression for the tunneling current! 31] 

/l^2 = ^ J^[/(^m,/il)-/(K,M2)]<5(^„^-£;„)|M„„|2, (30) 

m,n 

where E^ and En are the energy levels of eigenstates m and n on each side with the (spin 
excluding) wavefunctions ?/^m and respectively In a first order perturbation theory the 
transition matrix element was shown to be Il32ll 



I ^mn I 



dS ■ [V-mVV; - V'nVV'r, 

So 



2m ^ 

where the surface integral is carried out over a median surface So- 

If we define dimensionless wavefunctions ip, instead of normalized ones ip, by 



(31) 



fm{r) = V^iil^m{r), (32) 

then we can further write 



^ - ' / dS • [ifm^lfl - flV(pr> 

JSa 



'So 

and define a volume-independent transfer matrix element as 



(33) 



\2m J 



2 . 2 
/ dS • [(/7m V(/?* - v?* V(^^ 



(34) 



which appears in eq. [29] and fulfills the correct dimensionality of the integral. 

We apply the above formalism for calculating the tunneling current for the three tunneling 
cases in our system: (i) between the STM tip and the silicon sample; (ii) between the STM tip 
and an isolated DB; and (iii) between an isolated DB and the silicon sample. Below we adapt 
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the general formula above to each of these cases in order to derive the simpler forms given in 
the text. 

For the first case, we assume that the tip wavefunctions are given as in the Tersoff- 
Hamann approach by eq. [6] and its LDOS is constant in energy. Also, we assume that the 
sample LDOS is constant in space and has the customary square root dependence on energy 
for the bands in semiconductors 



9Si{E) = ^ml%^E-EcBM (35) 



for the conduction band levels (above CBM) and 

9Si{E) = ^^j^ rn^p%V Eybm - E (36) 

for the valence band levels (below VBM), with gsi{E) being zero in the band gap (mn,efr and 
rrip eflf are the effective masses for electron and hole in silicon). 

For the tunneling between a DB and an electrode (tip or sample), we take advantage of 
the fact that the LDOS of the DB is a (5-function in energy, namely g-DB{E) = 5{E — -E'db)/^^db, 



and therefore the integral over energy in eqj29] is reduced to a single-energy term for -Edb- 
This yields a simplified formula 

Aire 1 ~ 

-^tip-DB = — ^7;; Ptip(-E'DB)[/tip(-E'DB) — /DB]l^mn(-EDB)| , (37) 

n i^DB 

where JIdb is the volume of the DB, defined in our case as the region in space where the 
DB charge density is greater than 1% of its maximum. A completely similar equation can 
be written for the tunneling between DB and Si just by replacing the subscript "tip" above 
with "Si". In this latter case, Mmn depends on the tunneling barrier between the DB and 
CB states lying at the same energy, which is given by the form of the band bending in the 
Si region surrounding the DB. Using the semiclassical approximation, we assume the tails of 
the wavefunctions for states in the conduction band in the barrier region scale as 9506(2;) = 
exp(— CcB-z), with Ccb = \/2mU\^/h and C7bb the average tunneling barrier height between DB 
and CBM. 

For all cases described above, the tunneling matrix elements were calculated numerically 
by assuming the corresponding wavefunctions given in the text. 

B. CALCULATING BAND BENDING IN THE VICINITY OF A DB DURING STM IMAGING 



Tip-induced band bending was calculated within the semiclassical approximation using 
finite element methods (FEM). We assume that the STM tip has an axially symmetric shape 
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ending in a hemispherical surface of radius Rup = 20 nm. On top of this surface we assume 
the existence of an atomic-sized protrusion of radius Rq, which coincides with the spherical 
tip radius assumed in the Tersoff-Hamann model for the tip wavefunction. 
For an STM system, the total potential drop between the tip and the sample 

(Vtip - Kam) = 14ias + "^tip - X " (^CBM - Ep) (38) 

where ^tip and x are the tip workfunction and electron affinity of Si. Another component of 
Ubb can arise in the vicinity of a charged DB from the field created by the localized electron. 

The calculation involves solving a system of non-linear equations consisting of the Poisson 
and Schroedinger (PS) equations for the electrostatic potential U and the charge densities 
(n,p) of electrons and holes, respectively, in the semiconductor. Poisson equation reads 

V^U{r) = - — [N++p{r)-NX-nir)], (39) 

where e is the elementary charge, esi is the electric permitivity of the semiconductor, the 
ionized donor concentration, and the ionized acceptor concentration. The same U that 
satisfies the above equation, must also satisfy the Schroedinger equation 

- 7^V2^(r) + U{r)i>{r) = E^{v) (40) 

where, in our case, iJj are extended states of the crystal belonging to the CB or VB bands, and 
their corresponding charge distribution determines n and p. 

In the semiclassical approximation, the extended states ijj entering the above equation 
are affected by a local potential just by a rigid shift in their energies equal to —eU{r). As a 
consequence, the local charge densities are shown to be proportional to the (l/2)-Fermi-Dirac 
integrals, F1/2 



n(r) = iVeFi/2((^F - + eU{v))/kBT) 



(41) 



p(r) = iV.Fi/2((i?buii^ - eU{T) - Ep)/kBT) 

where E-p is the Fermi level (chemical potential) of the crystal. The above PS equations were 
solved self-consistently for the tip-induced band bending (no charge present on a DB) using 



an iterative FEM scheme in which the solution to eq. 41 was input as a source term in the 



RHS of Poisson eq. |39) To ensure convergence of the iterations the Anderson mixing scheme 
was used. 

In order to calculate the distortions in the bands caused by the field of a DB , we performed 
a second FEM calculation in the presence of the charged DB. For practical reasons, we did 
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not employ FEM on a three-dimensional geometry (for the case when the DB is laterally 
displaced from the central axis of the tip). Such a calculation is not practically feasible for 
the following reasons: (i) the size of the finite element domain must extend very deep into 
the semiconductor bulk in order to ensure that boundary condition there (zero electrostatic 
potential) does not introduce non-physical artifacts in the results; (ii) a very fine mesh is 
needed for the tip vicinity where the field is very strong in contrast to the bulk semiconductor 
where the field is relatively weak; (iii) many PS iterations are required to reach convergence, 
depending on the biasing conditions and doping level in the semiconductor. 

We found a way to approximate the full solution of the PS system in the presence of a 
charged DB by taking advantage of the fact that the STM tip has a radius much greater than 
the size of the DB, and therefore it can be seen as locally flat in the vicinity of the DB. This 
implies that the potential due to the DB alone, C/DB-(r), will be approximately the same (as 
measured from the center of the DB) whether the DB is directly under the tip apex, or slightly 
displaced laterally, by a few nm. This allows us to extract t/oB- (r) as the difference between 
two different FEM solutions for two axially symmetric systems: (i) STM tip alone, and (ii) 
STM tip with a DB directly under the apex. Then, for the case when the DB is not on the 
central axis of the tip but close to it, the solution is well approximated by adding the potential 
C/DB-(r — tdb) to the tip-induced potential. More details of the finite element calculation will 
be published in an upcoming paper [25]. 

The calculation of the charge carrier densities at the surface including charge quantization 
in the surface-perpendicular direction was also performed. We numerically solved the Pois- 
son and Schroedinger equations self-consistently for a planar metal-vacuum-semiconductor 
system in a one-dimensional case, assuming translational symmetry in xy-plane, in which 
case the Poisson and Schroedinger equations need only to be solved in the z-direction. This 
approximation is justified by the fact that the surface-parallel length-scale of the sample in 
which tip induced band bending is manifest is on the order of hundreds of nanometers, and 
much greater than the surface-normal length-scale over which the hole density is signifi- 
cant (tens of nanometers). In other words, the confinement of the holes in the potential well 
created by the biases STM probe in the surface-parallel direction is much weaker, than the 
confinement in the surface-normal direction. Charge quantization was also estimated using a 
ID Poisson solutionpMl corresponding to the limit of a very large tip radius (or flat tip); those 
results were consistent with ours. 
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C. STM VICINITY CAPTURE OF TUNNELING ELECTRONS BY A DB 

The rate of capture for a deep-level state is proportional to the velocity of the electron being 
captured[34], Cn = cr^VnU. In equilibrium, at finite temperature, electrons sit at the bottom of 
the conduction band, and their energy distribution is well described by Maxwell-Boltzmann 
statistics, as long as Ecbm — Msam > ^b^. This allows the capture rate to be described by 
the familiar equation Cn = a^vthn, expressed in terms of the thermal velocity. The thermal 
velocity is the average group velocity for all electrons in the conduction band. Assuming equi- 
librium statistics and a parabolic conduction band minimum, the thermal velocity is given by 
vth = x/SfeBT/vrm* ^ lO'^ cm/s. 

However, in calculating the excess capture rate, c^^, due to electrons injected from the 
tip into the conduction band, we no longer have recourse to Maxwell-Boltzmann statistics. 
We therefore use the group velocities directly. For simplicity, we retain the assumption of a 
parabolic band, so that -Ek = h^k'^ /2m. In addition to simplifying calculations, this ensures 
that group velocities are parallel to their associated wavevectors. 

The injected electrons can tunnel into eigenstates in the conduction band with energies 
below /Utip, however we expect that tunneling will occur most readily into states with high 
fcz (surface normal) and low K (surface-parallel) values. Thus their velocities are not evenly 
distributed across the range of polar angles, 6. We solve for the angular distribution of tip- 
injected electrons by considering the matrix element, |Mmn(k)p, for tunneling from tip states 



to sample states with wavevector, k, described by equation 34 The angular distribution of 



wavevectors, and hence also of velocities, of the injected electrons is then given by 

/■V^m(Mtip--BcBM)/^ 

D{e)=M j k^\Mmn{^)\ dk (42) 

where 







M = — . (43) 

Ultrafast pump-probe reflectivity measurements of the Si(OOl) surface place the momen- 
tum relaxation time of free carriers at 32 fs ll35]| . Electrons traveling with speed vth (much 
slower than the average velocity for injected electrons) will travel a distance of roughly 6 nm 
in this time. We therefore make the approximation that injected electrons retain their initial 
group velocities over the distance scales relevant to the present problem. The local excess 
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electron density at each point in the silicon sample, due to the injected current from the STM 
tip when the DB is neutral, I^p°si' then given by 

T-DB" 

nS™(r) = i^(e)_^ (44) 

where r = {r,9,(f)) being the distance vector from the STM tip apex, with the polar angle 9 
being measured from the tip axis pointing toward the silicon crystal. 

Finally we account for the occupation of the DB with one electron, 1 — /^g, and write the 
capture current by the DB level as 

Cn^™ = (1 - /5b) / dr IV^DBWr D{e)lE^lsi^„ (45) 

where we assume that the capture cross section for an infinitesimal volume is |^db|^ fnC^r- 
The capture rate in s~^ is obtained by dividing ^stm ^j^g elementary charge and 
assuming a neutral DB, i.e. /^g = 0. 
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